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Abstract 

This work proposes an explanation of the Pioneer anomaly, the un- 
modelled and as yet unexplained blueshift detected in the microwave 
signal of the Pioneer 10 and other spaceships by Anderson et aim 1998. 
What they observed is similar to the effect that would have either (i) 
an anomalous acceleration ap of the ship towards the Sun or (ii) an 
acceleration of the clocks at = ap/c. The second alternative is inves- 
tigated here, with a phenomenological model in which the anomaly 
is an effect of the background gravitational potential ^f(t) that per- 
vades all the universe and is increasing because of the expansion. It 
is shown that 2a t = d^/dt = d 2 r c i oc k s /di 2 , evaluated at present time 
to j where t and r c i oc k s are the coordinate time and the time measured 
by the atomic clocks, respectively. The result of a simple estimate 
gives the value a t ^ 1.8 x 10~ 18 s _1 , while Anderson et al suggested 
at = (2.9 ± 0.4) x 10~ 18 s _1 on the basis of their observations. The 
calculation are performed near the Newtonian limit but in the frame 
of general relativity. 

KEY WORDS: Pioneer anomaly, Pioneer acceleration, acceleration of 
the clocks, background gravitational potential 
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1 INTRODUCTION 

1.1 The anomaly. Anderson et al reported in 1998 the observation of 
an unmodelled Doppler blueshift in the microwave signals from the Pioneer 
10/11, Galileo and Ulysses spacecrafts that increases linearly in time 1 - 1 ). They 
had been observing it since more than twenty years. Obviously, its simplest 
interpretation is that the ships were not following the predicted orbits, as if 
our star pulled a bit too much from them with a force independent of the 
distance. The corresponding anomalous acceleration, directed towards the 
Sun and constant, would have the value^ 

a P = (8.74 ± 1.33) x l(T 10 m/s 2 . (1) 

Intriguingly enough, the effect does not show up in the planets. In their 
first paper they said: "it is interesting to speculate on the possibility that 
the origin of the anomalous signal is new physics" ^ 2 ) , and later "The ve- 
racity of the signal is now undisputed, although the source of the anomaly, 
some systematic or some not understood physics, is subject to debate" ®. 
For an interesting argument showing that it may not necessarily be due to 
systematics, see reference 4. The effect is still unexplained 1 - 2 ) . 

Anderson et al say that their data show*- 1 -* "a steady frequency drift of 
about —6 x 1CT 9 Hz/s, or 1.5 Hz over 8 yr. This equates a clock acceleration, 
—a t , of —2.8 x 1CT 18 s/s 2 ", what would mean that the frequencies would drift 
as 

v = v [l + 2at(t-t )}, (2) 

to being here the initial time of the observations (with the best value for ap 
a t = (2.9 ±0.4) x 10~ 18 s _1 ). The relation with the Pioneer acceleration 
is ap = a t c. This is important since they say that the drift in the Doppler 
residuals cannot be removed without either an acceleration of the ship ap or 
the inclusion of a "clock acceleration" a t . Such acceleration a t would imply 
that all the clocks would be changing with a constant acceleration or, in other 
words, that there would be a nonuniformity of time. They found that the 
first alternative leads to problems with the equivalence principle and with 
the cartography of the solar system. They considered the second by means 
of several models in which the time is distorted phenomenologically (see 
reference 2, sections XI.D and XI.E). The best results were obtained with 
a model that adds a quadratic term to the definition of the International 
Atomic Time. However, they found some problems and concluded: "The 
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orbit determination process clearly prefers the constant acceleration model, 
ap, over the quadratic in time model." 

1.2 Purpose and assumptions of this work. This paper considers this 
second alternative. It shows that, because of the expansion, the background 
gravitational potential that pervades all the universe produces an accelera- 
tion of the cosmological proper time with respect to the coordinate time. In 
its turn, this implies an acceleration a t of the atomic clocks. A simple esti- 
mation gives a value for a t that is close to that found by Anderson et al. The 
anomaly is thus in this model an effect of the expansion of the universe. It is 
assumed, for simplicity, that (i) all the matter and energy of the universe are 
uniformly distributed, (ii) the space sections t = constant are flat, and (iii) 
the near Newtonian approximation is adequate and meaningful. The time 
coordinate t is so chosen as to go in the Newtonian limit to the Newtonian 
time. The model here presented is the relativistic version of a previous New- 
tonian one 1 - 5-7 ), see also reference 8. Two other models in which the anomaly 
is also due to the expansion are proposed in references 9 and 10. 

1.3 Two definitions of the light speed. It is important to know precisely 
which one of the several meanings of "speed of light" is used^ 11 **. Particularly, 
it must be reminded here that the light speed can be defined in general 
relativity in two different ways, (i) with respect to the cosmological proper 
time r, c* = d£/dr (= constant), and (ii) with respect to the coordinate time 
t, c = d£/dt = c(r, t), where d£, dr and dt are elements of spatial distance, 
proper time and coordinate time along a null geodesic, respectively (the 
first is the usual definition). These two speeds will be denoted as c* and c, 
respectively, and will be called "proper speed of lighf and "ordinary or non- 
proper speed of light. The derivative of c with respect to t at present time t , 
denoted as ae = c(t ) = dc(to)/dt, will be called non-proper acceleration of 
light or just acceleration of light if there is no risk of confusion. The first is 
a universal constant of nature (as it must happen in general relativity), the 
second is not but, quite on the contrary, it depends generally on space and 
time c = c(r, t). The duality between c* and c reflects the relation between 
the proper time and the coordinate time. 

The element of interval can be written* 12-15 ) (assuming for simplicity that 
9oi = 0) 

ds 2 = c* 2 dr 2 - df , (3) 
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with dr = i/^oodt and d£ 2 = gijdx l dx^ so that c* is constant and c is equal 
to 

c = c (r,t) = c*^/5oo. (4) 

Near the Newtonian limit, (700 — 1 + 2$/c 2 , at first order, $ being the 
gravitational potential, so that c = c(r, t) = c*[l + $(r, t)/c 2 ]. Taking a non- 
zero origin for the potential at a reference laboratory R, this can be written, 
at first order, as 

c = c ( r , t) = c [l + $(r, t)/c 2 (r, £) - $ fl /c 2 ] , (5) 

where Cq and are the values of c and the potential at that laboratory. 

In this paper, we are interested in the effect of the background gravita- 
tional potential that pervades all the universe and is due to all the existing 
matter and energy. Assuming a uniform distribution of matter and energy, it 
is clear that is time depending but space independent. It will be denoted as 
$av(£) ("av" stands for average since the mass-energy density is averaged). 
The same must happen therefore to c = c(t). Instead of © one would have 
then 

C = C(t) = C [l + $ av (i)/ C 2 (t) - $ av (t )/Co] , (6) 

with Co = c(t Q ), to being a reference time that will be the present time or 
age of the universe, in general. Because of the expansion, it turns out that 
$ av (t)/c 2 (t) is an increasing function, as will be shown in section 4 where 
its derivative with respect to t will be calculated. Consequently, c(t) is also 
increasing. 

If the universe would contain only matter, be it ordinary or dark, the 
background potential $ av would be negative so that 

C(t) < C [l - $av(t )/^] • (7) 

As will be seen in section 3, this implies c(t) < c*, as it could be expected. 

However, the effect of the cosmological constant or of dark energy changes 
dramatically this question. The cosmological model used here is the standard 
with 27 % of matter and 73 % of dark energy. For the latter we take either 
a cosmological constant or the quantum vacuum, but the conclusion would 
be the same for any kind of dark energy with an equation of state implying 
repulsion. The potential $ av (£) is the addition of the two effects of matter 
and dark energy. 
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What is important here is that both the cosmological constant and the 
quantum vacuum produce a positive potential. Furthermore, it turns out that 
$ av (t) is an increasing positive function after a certain time, in particular now 
(this will be shown in section 4). Equation (JHJ) implies then that c increases as 
far as the universe expands and, in particular, it can be larger than the proper 
speed of light c*. Although this maybe seem strange and contrary to current 
wisdom, it must not be a matter of concern since, as it must be emphasized, 
the proper light speed of light c* is in fact a universal constant in this model. 
Moreover, c < c* if there is only matter, let it be ordinary or dark. It is only 
the dark energy, which, in addition to accelerate the universe, can make that 
c(t) could be larger than the proper speed of light c*. This intriguing result 
will be considered in section 3. To understand this property, one must keep 
in mind that almost all our intuitions in general relativity were generated in 
the study of the gravitation of matter (i. e. without dark energy). 

There is moreover a functional relation between the two times r = r(t), 
such that t must accelerate with respect to t. Because of the freedom to 
choose the coordinates in general relativity, this last statement must be qual- 
ified: the time coordinate t is defined here so that in the Newtonian limit it 
goes over the Newtonian time. All this means that 

cm = — = c — and — - — = c —rr . (8) 

w dr dt dt dt dt 2 y ' 

Equation (JHJ) states that the non-proper light speed c(t) must increase if the 
proper time r accelerates with respect to the coordinate time t, its time deriva- 
tive being equal to the proper light speed c* times the second derivative of r 
with respect to t. 

The expression "light speed" means usually the proper speed c* that, 
being a universal constant, is of the utmost importance. However, the non- 
proper light speed c is also used in some important cases. For instance, in 
the study of the bending of a light ray grazing the Sun surface. Let M and R 
be the mass and radius of the Sun. The interval around any star is given by 
the Schwarzschild metric, what implies that c = c(r) = c^l — r]R/r), with 
Coo = c(oo) and r\ = GM/c^R ~ 2.1 x 10 -6 . Einstein gave two formulae for 
this effect. The first (1907) is based only in the equivalence principle and 
gives (f) = 2rj = 0.875", just one half of the observed effect. The second (1916), 
in the frame of general relativity, gives the complete result = 4rj = 1.75". 
The first one can be obtained simply by considering the propagation of a 
wave light with the previous value of the non-proper light speed, in other 
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words as the solution of the variational problem 




2 1 + r)R/r 



d£ = 0, 



(9) 



Ji 



c, 



•oc 



where d£ = dx 2 +dy 2 +dz 2 is the Euclidean line element, i.e. as a consequence 
of the application of the Fermat principle to the non-proper light speed c. 
The complete effect is obtained by taking instead the non-Euclidean spatial 
line element of the Schwarzschild geometry. In other words, the problem is 
solved by assuming that the light propagates through space with the non- 
proper speed c(r) = 0^(1 — i]R/r) (taking into account the Riemannian 
character of the spatial metric). 

These considerations can be summarized as follows: because of the back- 
ground gravitational potential of all the matter and energy in the universe, 
(i) the cosmological proper time r accelerates with respect to the coordinate 
time, and (ii) the non-proper speed of light c(t) increases so that an > 0. 

1.4 Plan of the paper. In section 2, the Maxwell equations in general 
relativity are reviewed, with emphasis on the ideas of permittivity and per- 
meability of geometrical origin and on their effect on the non-proper speed of 
light c(r, t) . It is shown that 2a t must be equal to the derivative with respect 
to t of the background potential over c 2 (t) and that the corresponding non- 
proper of light acceleration at implies a blueshift. In section 3, the reasons 
for the acceleration of the atomic clocks are considered, showing that it is 
equal to 2a t . In section 4, an estimation of a t is carried out, the agreement 
being good with the result found by Anderson et al on the basis of their 
observations^. The conclusions will be stated in section 5. 



2 THE SPEED OF LIGHT 

AND THE MAXWELL EQUATIONS 



In order to understand the behavior of the non-proper speed of light c(r, t) 
(different, don't forget, from the proper speed of light c*, which is a universal 
constant), let us consider the effect of a gravitational field on the Maxwell 
equations, in which the time derivatives are with respect to the coordinate 
time t. As was seen in section 1.3, the proper speed of light c* and the non- 
proper speed of light c are equal in absence of potential or, in other words, 
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the reason for their difference is the presence of matter and energy. It is clear 
that near the special relativity or the Newtonian limits, the electromagnetic 
fields obey the classical wave equation with the velocity c(r, t) (remember 
that flat space sections are assumed). That local value of c must be used at 
any spacetime point in the wave equation. Let us take a reference laboratory 
where c = cq. The following discussion is based on the well known textbook 
The Classical Theory of Fields by Landau and Lifshitz (reference 16, section 
90). 

The electromagnetic tensor is defined in general relativity by means of a 
vector field such that F^ v = A u .^ — A^ u = d tl A u — d u A fl . The electromagnetic 
vectors E, D and antisymmetric tensors B^, H^ are defined as follows 

Ei = F Q i , Bij = Fij , D l = —y^gooF 01 , H lJ = ^/gooF^ , 

the vectors B, H being the dual to the three-tensors B^ and Hy, i. e. 
B i = _ e i J' fc J B jfc /(2 v /7), Hi = -^ye ijk H^ k /2, where 7 = detfry), 7^ = -g^ 
being the three-dimensional metric tensor (assuming for simplicity g^ = 0). 
It follows that^ 16 ) (eo = 1, /io = 1 in this argument) 

D = E/ V /^, B = H/v^. (10) 

If the space is empty, i. e. without free charges or currents, the Maxwell 
equations can be written as 

V-B = 0, VxE = -i^(V7B). (11) 

V-D = 0, VxH=l9 1 ( v ^D). (12) 

In a static situation, these four equations have exactly the same form as in 
special relativity, since the factors ^7 cancel. However, eq. (fTUj) implies that 
the relative permittivity e r and permeability /i r of empty space are differ- 
ent from 1, their common value being e r = /i r = (goo)~ 1//2 - This is due to 
the geometry of spacetime. If the potential depends on time and near the 
Newtonian limit one has g 00 = 1 + <3>(r, £)/c 2 (r, t) — $r/cq, where $r is the 
potential at a reference laboratory where c = Co, so that the empty space is 
like an inhomogeneous optical medium with 



e r (r, t) = Mr(r, t) = 1 - [$(r, t)/c 2 (r, t) - $ R /c 2 ] . (13) 
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Since c = co/y/e r fj, T , the non-proper speed of light is given as 

c(r,t)=Co{l + $(r,t)/c 2 (r,t)-$ R / C 2} , (14) 

$r being a reference potential, at present time in a terrestrial laboratory 
R, where the observed light speed is Cq. As a historical comment, this last 
equation was first obtained in 1911 by Einstein himself, as a first order ap- 
proximation in the static case, in a paper entitled "On the influence of grav- 
itation on the propagation of light" ^ (in 1907 he had already shown that 
c depends on $ as a consequence of the equivalence principle. Note that it 
is still valid in general relativity, at first order.) After a discussion on the 
synchronization of clocks, he concludes there "if we call the velocity of light 
at the origin of coordinates Cq, where we take $ = 0, then the velocity of 
light at a place with gravitational potential $ will be given as 

c = c (1 + $/c 2 ) ." (15) 

Einstein had not as yet introduced the general relativistic idea of proper time 
and used only the coordinate time to define the speed of light*- 18-19 -*. Note 
that (fH|) reduces to (fTKj) in the static case if the potential at the reference 
laboratory vanishes and that the consideration of the Maxwell equations 
confirms eq. (0), obtained in section 1.3 from the expression of the interval. 

In this work, we are interested in the case of a potential depending only 
on time. Let $ av (£) be the background potential of all the matter and en- 
ergy and ^(t) = $ av (t)/c 2 (t) its dimensionless expression (the subindex "av" 
being omitted in \1/ in order to simplify the notation). The permittivity and 
permeability (fT3*|) take then the form near present time to (>e. the age of the 
universe) 

6r (t) = ^(t) = 1 - - . (16) 

The non-proper speed of light is then 

c(t) = co [1 + *(*)- *(*„)] , (17) 

where Co = c(t ). This can be written as 

c(t) = c [l + 2a t (t - t )] = c + a £ (t - to), (18) 

the quantity a t and the non-proper acceleration of light ae = c(to) being 

2ot = *(to) , at = 2a t c = 2^(t )c . (19) 

As will be seen in the following, (i) is in fact 2ap, being therefore an 
adiabatic acceleration, and (ii) a t is what Anderson et al termed acceleration 
of the clocks. 
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2.1 The blueshift. It will be shown now that the non-proper acceleration 
at implies a blue shift of the light with respect to the coordinate time, at 
first order in ag. More precisely, it turns out that the frequency v of a 
monochromatic light wave with such an adiabatic acceleration ag increases, 
its derivative with respect the coordinate time t, z>, satisfying 

u/u = ai/c . (20) 

This means that an adiabatic non-proper acceleration of light has the same 
radio signature as a blue shift of the emitter, although a peculiar blue shift 
with no change of the wavelength (i.e. all the increase in velocity is used to 
increase the frequency). 

The derivative with respect to t of the background gravitational potential 
of all the universe \I/(t) is positive and of the order of the Hubble parameter 
Hq = 2.3 x 1CT 18 s , since the galaxies are separating (a calculation will be 
made in section 4). Equation (fTBj) tells then that the relative permittivity 
e r and permeability /x r of empty space are decreasing, their derivatives with 
respect to t being negative and also of order H , i. e. very small. This can 
be expressed by saying that the optical density of empty space is decreasing 
adiabatically. To study the propagation of the light in a medium with time 
depending permittivity and permeability, we must take the Maxwell equa- 
tions and deduce the wave equations for the electric field E and the magnetic 
intensity H, which are V 2 E - d t (M(eE)) = 0, V 2 H - d t (efl^H)) = 0, or, 
more explicitly, 

V 2 E-«9 2 E/c 2 (t) (/ t //i + 2e/ eo )^E/c 2 (t)-e/ i E/( eo/ u c 2 (t)) = 0, 

(21) 

V 2 H-«9 2 H/c 2 (t) (2/ t / / io + e/eo)ftH/c 2 (t)-e/ t H/( eo /ioc 2 (t)) = 0, 

(22) 

with c(t) = Co + a t (t — to), since at present time e r = 1, /i r = 1. Because 
e/e and /i/y^o are of order H = 2.3 x 1CT 18 s -1 , the third and the fourth 
terms in the LHS of (|2*T|) and ({22} can be neglected for frequencies ui 3> H , 
in other words for all practical purposes. We are left with two classical wave 
equations with time dependent light velocity c(t). 

V 2 E - <9 2 E/c 2 (t) = 0, V 2 H - <9 2 H/c 2 (t) = 0. (23) 

In order to find the behavior of a monochromatic light beam according to 
these two wave equations, we take for instance the first one and insert E = 
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E exp{— i[nz — (too + to{t — t )/2)(t — t )]}, where the frequency is the time 
derivative of the phase of E, i.e. uj + uj{t — t ). Neglecting the second time 
derivatives and working at first order in u (with uj{t — t ) <C Uq, Co <C u>q), 
substitution in (}2*T]) gives n 2 = [(lj + u(t — to)) 2 — iw\ /c 2 (^)- It follows that 
K = k + i( = ±(u /c(t))[l + uj{t — t )/^o](cos <p + % sin (p), with <p = —uj/2uj 2 , 
so that k = ±(c<j /c ) (1 + to(t — t )/u )/(l + a^it — to)/ c o) what implies 
k = ±ujq/c , Co/ujq = ai/c , as stated before. Equation (J2U|) has thus been 
proved. Also, ( = —Co/2ojqCq = g^/2cq. The wave amplitude decreases in the 
direction of propagation as e~ z l l with i = 2cq/g^, but as ai is of order H co, I 
is of order of 5,000 Mpc, or, in other words, this attenuation can be neglected. 
As is easy to show, to take k + kt for the wave vector leads to k = 0. These 
results are equally valid for the second equation in (|23|). Taking into account 
that a>i/c = 2a t and according to (fTo^ - (fT^|) . the frequencies drift as 

i/ = i/ [l + 2at(t-to)], (24) 

what shows that 2a t = ^(io) is the acceleration of clocks mentioned by 
Anderson et al. According to these arguments, its value is the derivative with 
respect to the coordinate time t of the background gravitational potential of 
all the universe. This will be further studied in next section. 

3 THE ACCELERATION OF THE ATOMIC 
CLOCKS 

3.1 The effect of the dark energy. Can c(t) be larger than c*? The 

observations are made by using atomic clocks, which measure proper time r 
not coordinate time t. However, up to now we have used mainly in this paper 
the coordinate time. This question is addressed in this section. Taking the 
t derivative of (fP7j) at time t , it is found that c(t ) = c(to)^(t ). The same 
argument can be applied to the expression for c(t) near any other fixed time 
t, what implies that c(t) = c(t) exp[^/(t) — ^(t)], an expression valid for Vt. 
In particular, taking t — to, one finds 

c(t) = Coe r*(*)-*(*°)]. (25) 

As can be seen, (II 7|) is the first order approximation to (|25|). The shape of 
the function c(t) defined by (|2*3|) does not change if the reference time (to or 
i) is changed, because c(ti)e - *^ 1 ^ = c(t 2 )e~*^ 2 ^. Since c* and c are equal if 



Acceleration of the clocks 



11 



^ = 0, it follows that c* = c e what does not depend on the particular 

value of to- In other words, the non-proper speed of light is 

c(t) = c*e*<*> . (26) 

We can precise now some ideas exposed in section 1.3. If there is only matter 
in the universe (i. e. without dark energy), then \I/(t) is necessarily negative, 
so that c(t) < c*. This is the usual situation that everybody has in mind. 
However, if there is dark energy, ^(t) can be positive, and in that case 
c(t) > c*. We find thus that the non-proper speed of light could be larger than 
the proper speed and constant of the nature c*. Even if this is a surprising 
effect of the dark matter, in addition to accelerate the expansion, it must be 
emphasized that it is compatible with the constancy of c* and with general 
relativity therefore. Indeed, it is an unexpected and new effect of the dark 
energy. 

3.2 On the time of the atomic clocks. The interval can be written then 

as 

so that 

dr = e* (t) dt and 



dt 2 



= ^(t )e* ( ' o) =2a t e* ( * o) >0, (28) 

to 



since \l/(t) is now an increasing function and if (to) > (the precise calculation 
is done in next section). Note that r is a well defined cosmological proper 
time and that it accelerates with respect to t, its second derivative being 
in general non nil. This means that if, sometime t\ in the past, an atomic 
clock was set to click at the same rate as the coordinate time, it would be 
advanced now with respect to the coordinate time. In particular, the two 
time intervals would be different now, since dr = e*(*°>df. If there is only 
matter, be it ordinary or dark, then ^(t ) < and dr < dt. On the other 
hand, if the dark energy exists, then ^(to) can be positive and then dr > dt. 
This second possibility is the one that actually happens: the background 
potential not only increases, but is positive also because of the effect of the 
quantum vacuum, or the dark energy, as is shown in next section. 

However, in the actual measurements these differences between dt and dr 
do not occur, since both times are based now on the international second, 
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defined with reference to the period of a transition of the cesium atom. This 
means that a small interval of the atomic clock time and of the coordinate 
time are equal, dr c i ocks = dt, at precisely t = t . Therefore, the time of the 
atomic clocks has been renormalized, indeed, multiplying the cosmological 
proper time by the constant scale factor e - *^ ^. This argument shows that 
the time really measured by the atomic clocks is 



dr, 



clocks 



,[*(t)-*(to)] df so that 



d 2 



^clocks 



dt 2 



4>{t Z 



2a t . 



(29) 



to 



This explains the real meaning of a t : it is one half the second derivative 
with respect to the coordinate time t of the time of the atomic clocks, if they 
have been renormalized by a multiplicative factor to tick now, just at time t , 
at the same rate as the coordinate time (in order for both to use the same 
definition of second). Note that there is a difference between the meanings 
of the "acceleration of the clocks": here it is 2a t while for Anderson et al. 
it is a t . Of course, in the future the atomic clocks will advance over t, i. e. 
dr c i ocks > dt if t > t . 

Taking the first order approximations near to, the intervals of both times 
are related as dr c i oc ks = [1 + ^(t) — ^(t )]dt. This is equal to dr c i oc k s = 
[1 + #(* )(* - t„)]dt = [1 + 2a t {t - t )]dt (compare with ©). It follows 



^clocks — Tclocks(^o) — [t — to) + Ot(t — to)' 



(30) 



what reminds the quadratic in time model tried by Anderson et al to solve 
the anomaly. Probably, this explains why this was the best among the other 
models they used with phenomeno logically distorted time (see section 1.1, at 
the end). 

In order to refer to the time of the clocks, a t must be multiplied by the 
factor dt/dr c i oc ks what produces the change 



a t =>• a 7 



dt 



dr, 



a t [l - 2a t (t - t )] = a t 



(31) 



clocks 



neglecting the terms of second order in a t . Since a t = a T , at first order, eq. 
(2) can be written as 



v = u {l + 2a T [r docks - r clocks (t )]} 



(32) 



at first order, what states that the frequency must drift, according to the 
measurements made with atomic clocks. This gives a solution to the riddle. 
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Indeed the two derivatives of the background potential \I/ with respect to 
the times t and r c i oc k s are equal because the intervals of the two times are 
equal at precisely the time t (or r c i oc k s (to))- Consequently, the acceleration 
of the clocks can be calculated as one half the derivative of ^ with respect 
to any one of the two times (at first order). Therefore, this model seems 
to offer an attractive solution to the anomaly. But in order to test it, its 
prediction for the inverse time a t must be calculated and compared with the 
value proposed by Anderson et al. This is done in next section. 

4 ESTIMATION OF THE ACCELERATION 
OF THE CLOCKS 

Unfortunately, a rigorous calculation of a t that would take into account all 
the eventual effects, is not easy. However, a simple crude estimate of its value 
can indeed be made, as is done in this section, taking (|T7j) and (jTUj) as starting 
point. Although involving approximations and simplifications, the result is 
meaningful since it shows the main ideas of the model, giving a convincing 
representation of the phenomenon. In order to do that, the shape of the 
function ^(t) near to must be determined first to calculate its time derivative 
eq. The potential of all the universe at the terrestrial laboratory R can 

be written, with good approximation, as $ a n = $i oc (-R) + $av(*)- The first 
term $i oc (i?) is the part due to the local inhomogeneities, i. e. the nearby 
bodies (the Solar System and the Milky Way). It is constant in time since 
these objects are not expanding. The second $ av (t) is the space averaged 
potential due to all the mass and energy in the universe (except for the 
nearby bodies), assuming that they are uniformly distributed. Contrary to 
the first, it depends on time because of the expansion. The former has a non 
vanishing gradient but is small, the latter is space independent, but time 
dependent and much larger. The value of $i oc /co at the laboratory R is the 
sum of the effects of the Earth, the Sun and the Milky Way, which are about 
—7 x 10~ 10 , — 10~ 8 and —6 x 10~ 7 , respectively, certainly with much smaller 
absolute values than $ av (i ), of the order of — 1CT 1 as will be seen below. 

Since the background gravitational potential of all the universe $ av (t) is 
increasing because of the expansion (the galaxies are separating and their 
interaction potential increasing), eqs. (jl9j) imply that a t and ag are both 
positive, as will be seen in the following. In this sense, there is a non-proper 
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acceleration of light ag = c, see section 1.3 (as it may be convenient to stress 
again, all this is compatible with the constancy of the proper speed of light 



Let fijvi = 0.27, = 0.73 be the corresponding present time relative 
densities of matter (ordinary plus dark) and dark energy corresponding to 
the cosmological constant A. We take a universe with k = and Hubble 
parameter H = 71 km • s" 1 ■ Mpc -1 = 2.3 x 1CT 18 s -1 . In order to deter- 
mine the average potential $ av (t), let $o(^o) be the gravitational potential 
produced by the critical density of mass distributed up to the present ra- 
dius of the visible universe -Ru(^o) = c o/H = 4,200 Mpc. One has then 
$ (t )/ c 2 = - J o co/Ho Cq 2 Gp cr 47rrdr = 2ixGp cv /Hl ~ -0.78. It must be em- 
phasized that, although this value of the potential might seem to be too large 
for this approximation to apply, there is no problem in fact since it is space 
independent and its time derivative is extremely small. It will be absorbed 
in the redefinition of time, its effect being only to accelerate adiabatically 
the proper time with respect to the coordinate time, i.e. to accelerate the 
clocks, precisely what we want to investigate. 

Because the radius of the universe is changing, the potential must be 
multiplied by the factor [i?u(t)/i?u(to)] 2 , with R\j(t ) = c . It turns out then 



The present time value of the background potential is then & av (t ) = 
$ (to)(^M — 2^a) ~ 0.92 > (it is positive because of the contribution of 
the cosmological constant). Because of the expansion of the universe, the 
gravitational potentials due to matter and dark energy equivalent to the cos- 
mological constant vary in time as the inverse of the scale factor R(t) and as 
its square R 2 (t), respectively (with R(t) = (fi M /fi A ) 1/3 sinh 2/3 [(3A) 1/2 t/2] 
for this model universe). This implies that the average background gravita- 
tional potential can be expressed as 



(remember: \l/(t) = $ av (^)/c 2 (t))- Note that, as announced in section 1.3 and 
since $o(^o) < 0, the term in fl^ overcomes necessarily the one in fl M , after a 
certain time, the potential being positive afterwards, as far as the expansion 
goes on. The non-proper speed of light c(t) is then larger than c*. As is easy 
to see, $ av (t) — > at time zero, because R\j(t) ~ t 2 in that limit. After a bit 



c* = di/dr.) 



that %{t)/cl = 27rGp cr R{ ] (t)/cl 




(33) 



Acceleration of the clocks 
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of simple algebra, the inverse time a t = ^ r (to)/2 (jTHj) can be expressed as 

Ho (1 - 9O A )$ /cg 

at 2 i + 2(i-3ft A )VqT 

Introducing in this equation the values of Qa and $o> the acceleration of the 
clocks a t turns out to be 

a t ~ 0.8 # - 1-8 x KT 18 s" 1 . (35) 

This is to be compared with the value suggested by Anderson et al on 
the basis of their data a t = (2.9 ± 0.3) x 10 -18 s _1 . Taking into account the 
simplicity of the calculation and the approximations involved, the prediction 
of this model can be considered to be acceptable. This is encouraging. 



4.1 An intuitive and phenomenological understanding of the phe- 
nomenon. Let a photon with frequency hv$ travel across the space at time 
r . Its energy is hu . If the time runs from r until r, when the gravitational 
potential is \I/(t), it will pick up potential energy, its total energy being 
then hv{r) = /iz/ [l + ^(t) — ^(to)]. It will be seen as having a frequency 
u{t) = vq [1 + ^(r) — ^(ro)]. The derivative of the frequency with respect 
to the time of the atomic clocks will be 2a T = d\l//dr c iocks > 0. In other 
words, it must be expected that the frequency of the photons will increase 
adiabatically, because of the expansion of the universe. This means a small 
blueshift, just what was observed. Since c* is constant, the wavelengths must 
decrease accordingly. That is the Pioneer anomaly. 



5 SUMMARY AND CONCLUSIONS 

The model here presented suggests an explanation of the Pioneer anomaly 
that is simple and based in standard physical ideas: the acceleration ap is not 
related to any anomalous or unmodelled motion of the spaceships. Instead, it 
is an effect of the increasing background gravitational potential that pervades 
the universe and produces an acceleration 2a t of the time of the clocks r c i oc k s 
with respect to the coordinate time t, i. e. a t = d 2 r c i oc k s /dt 2 | <o /2, to being the 
present time (remember that the time coordinate was so chosen as to go to 
the Newtonian time in the Newtonian limit). This acceleration of the clocks 
2a t is also equal, in this model, to the time derivative of the background 
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gravitational potential ^(i), i. e. a t = d$/dt\ t /2 = d^//dr c i O cks| t0 /2 (note 
that dr c i oc k/dt = 1, at present time t but not before or after). The anomaly 
would be thus an interesting case of the dynamics of time 1 - 20 ' 21 ). A further 
comment: as it might be worth to point out, this model is similar to the 
explanation by Mach of the origin of the inertia. 

According to this model, the anomaly is a manifestation of the expansion 
of the universe, which causes the increase of the background potential ^(t). 
This increase, in its turn, accelerates the cosmological time and causes the 
acceleration 2a t . A simple estimation gives a good agreement with the value 
proposed, on the basis of their observations, by Anderson et al, the discoverers 
of the anomaly (see eq. (J3SJ).) 

The model complies with the principles of general relativity, particularly 
because the proper speed of light (i. e. c* = d£/dr, r being the cosmological 
proper time) is a universal constant. However, what is here called the non- 
proper speed of light (i. e. c(t) = d£/dt, t being the coordinate time) is not 
constant. This is standard. If there is only matter, be it ordinary or dark, 
then c(t) < c*. On the other hand, it happens that c(t) can accelerate until 
being larger than the proper light speed c*, if and only if there is cosmological 
constant (or any other form of dark energy that implies repulsion). This is 
unusual and might seem perplexing at first sight, but it does not imply any 
contradiction as far as c* is still constant. It would be just another unexpected 
effect of the dark energy, in addition to accelerate the universe. 

The Pioneer anomaly poses a most intriguing riddle for physics and very 
complex and difficult problems for metrology. Taking this into account, the 
main conclusion of this paper is that the ideas here presented should be con- 
sidered by the experts who know the details of the motion of the spacecrafts 
and of the metrological procedures involved in the observation. 
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